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ABSTRACT 

The Seyfert 1 galaxy Ark 564 was observed with Chandra high energy trans- 
mission gratings for 250 ks. We present the high resolution X-ray spectrum that 
shows several associated absorption lines. The photoionization model requires 
two warm absorbers with two different ionization states {logU = 0.39 ±0.03 and 
logU = —0.99 ±0.13), both with moderate outflow velocities (~100 km s~^) and 
relatively low line of sight column densities {logNn = 20.94 and 20.11 cm~^). 
The high ionization phase produces absorption lines of O Vll , O Vlll , Ne ix , 
Ne X , Mg XI , Fe xvil and Fe xvill while the low ionization phase produces lines 
at lower energies (O VI & O vil ). The pressure-temperature equilibrium curve for 
the Ark 564 absorber does not have the typical "S" shape, even if the metallicity 
is super-solar; as a result the two warm-absorber phases do not appear to be in 
pressure balance. This suggests that the continuum incident on the absorbing gas 
is perhaps different from the observed continuum. We also estimated the mass 
outflow rate and the associated kinetic energy and find it to be at most 0.006% 
of the bolometric luminosity of Ark 564. Thus it is highly unlikely that these 
outflows provide significant feedback required by the galaxy formation models. 
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Introduction 



Outflows are ubiquitous in AGNs, manifested by high-ionization absorption lines in X- 



rays and UV (lReynoldslll997l . ICrenshaw et a 



related to the accretion process (e.g., iProga 



20031 . and references therein) and are perhaps 



20071 ). Understanding outflows is therefore as 



important as understanding accretion itself. The X-ray absorbers, commonly known as warm 
absorbers (WAs), have typical ionization parameter log ^ of — 2 erg s~^, a column density 
of Nh = 10^° — 10^^ cm~^ and an outflow velocity of 100 — 1000 km s~^, produced by warm 
ionized gas (T ~ 10^ — 10^ K; iKrongold et al.l 120031 ). The WAs are usually detected in the 
0.3 — 2 keV soft X-ray band by absorption lines of Oxygen (Ovill , Ovil , Ovi ), Iron 
(Fe Vll-xil and Fexvil-xxil ), and other highly ionized elements. Transitions by Civ , N v , 
and O VI are observed in both the X-ray and UV spectra of these sources with similar outflow 
velocities, suggestin g a connection bet^yeen the narrow absorption line systems in th e UV 
and the X-rays WA jMathur et al.lll994j : Il995l . iKaspi et allbood iKrongold et~aDl2003h . 



Several phenomenological models have tried to explain AGN warm absorber spectra 
showing multiple velocity components of multiple lines, and after years of effort a consensus 
is growing. In majority of the cases, if not all, physical properties and kinematics of the 
absorber are well determined and it can be described by at least two discrete ionization 
components (jPetmers et al.ll201ll . iHolczer fc Behaiil2012l ). These components are consistent 
with the same outflow velocity and appear to be in pressure equilib r ium, and so likely 
emerge from a mul t iphase wind (e.g., Krongold et al.l l2003l : l2005l : 120071 . iNetzer et al.l 12003 



Cardaci et al.l l2009l . lAndrade- Velazquez et al.ll2010[ ). The low-ionization phase (LIP) of the 
wind produces UV and X-ray absorption lines, but the high-ionization phase (HIP) is seen 
only in X-rays. 

Despite ubiquitous detection and successful modeling of WA spectra by multiple ab- 
sorbing components, very little is still known about their geometry and dynamical strength. 
Where do these outflows originate? Proposed locations span a wide range, of a factor of 10^ 
in radial distance from the central ionizing source: the accretion disk (as s uggested by the 
accretion disk wind models; 



20051 ) , the obscuring torus (iDorodnitsvn et al. 



Proga fc Kallman 



and to the narrow line region (jBehar et al. 



200 



2008! 



2003 



), the broad line region (iKraemer et a 



Crenshaw et al 



20o3. In principle, these 



Krolik fc KrisjboOlllBlustin et al.lboosi ) 



outflows could potentially provide a common form of AGN feedback required by theoretical 



models of AGN-galaxy formation fISilk fc Reed Il998l . IWvithe fc Loebl l2003l . iFabian 1 12012 



^The ionization parameter ^ — Lion/nr'^, where Lion is the ionizing luminosity between 1 Ryd and 1000 
Ryd (1 Ryd=13.6 eV), n is the number density of the material and r is the distance of the gas from the 
central source. 
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and references therein), although estimating their mass outflow rate and the kinetic energy 
outflow rate depends critically on the location of WAs. Some recent studies of WAs found 
that the typical outflow velocity is a small fraction of the escape velocity and that WAs do 



not carry sufficient mass/energy/ momentum to be efficient agents of feedback (IMathur et al. 
20091 . iKrongold etaPboO?! : boiol ). 



With the goal of self consistent analysis and modeling of grating spectra of WAs, we 
present here the results of our analysis of Chandra archival data of Ark 564. The analy- 
sis of one Chandra obse rvation of this source done in 2000 (50 ks) has been published by 
Matsumoto et al.l (120041 ): here we present new data of the three 2008 Chandra observations 
(250 ks total). 



Ark 564 



Ark 564 is a bright, nearby, narrow-line Seyfert 1 (NLSl) galaxy, with z = 0.024684, 



V = 14.6 mag (de Vaucouleurs et a. 



( iTurner et al.ll200lUMatsumoto et a. 



wave bands (e.g. 



Turner et al.l 



2001 



1991), and L2-10 kev = (2.4 - 2.8) x 10^^ ergs s"^ 
2OO4J. and present work). It has been st udied across all 



Collier et al. 2001, Shemmer et al. 2001, Romano et al. 



20041 ) and shows large amplit ude flux variations on short time scales and a peculiar emissio n 
line-like feature n ear 1 ke V ( iLeighly et al.l Il999l . iTurner et al.l I2OOII IComastri et al.l I2OOII ) . 
Matsumoto et al.l analyzed the Chandra HETGS observation of Ark 564 (that carried out in 
2000) with an exposure time of 50.2 ks. They modeled the hard X-ray spectrum with a power 
law of photo-index of 2.56 ± 0.06 and fit the soft excess below 1.5 keV with a blackbody, of 
temperature 0.124 ±0.003 keV. They find some evidence for a two phase WA with ionization 
parameters log ^ ~ 1 and log ^ ~ 2 and column density of logN^ = 21 cm~^. They find 
that the 1 keV emission feature is not due to blends of several narrow emission lines and 



suggest that it could be an artifact of the warm absorber. iBrinkmann et al.l ( 120071 ) studied 
the spectral variability of the X-ray emission of the Ark 564 using the ^ 100 ks XMM- 
Newton observation and find that the "power law plus bremsstrahlung" model describes the 
spectrum well at all times, with fiux variations of both components. 



Papadakis et al.l (120071 1 analyzed the XMM-Newton EPIC data from Ark 564 2005 



observations. They found evidence for two phases of photoionized X-ray absorbing gas 
with ionization parameter logC, ~ 1 an d log ^ ^ 2 and column densities of A^^^^ ~ 2 and 
5 X 10^° cm~^, similar to the results of iMatsumoto et al.l ( 120041 ). They also detect an ab- 
sorption line at ~ 8.1 keV in the low resolution CCD spectra and assuming that this line 
corresponds to Fe xxvi Ka, they suggest the presence of highly ionized, absorbing material 
of Nfj > 10^^ cm~^ outfiowing with relativistic velocity of ~ 0.17c. 
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AnalyzinK both th e EPIC and RGS data from the same XMM-Newton 2005 observations, 
Dewangan et al.l (120071 ) find two warm absorber phases with ionization parameters log ^ ~ 2 
and logC, < 0.3 and column densities of Nh ~ 4 and 2 x 10^*^ cm~^ and outfiow velocities of 
300 and 1000 km s~^ respectively. 



Smith et al.l (120081 ) analyzed the combined XMM-Newton RGS spectrum of Ark 564 



obtained from 2000 to 2005. They found three separate phases of photoionized X-ray ab- 
sorbing gas, with ionization parameters log ^ = —0.86, 0.87, 2.56 and column densities of 
Nh = 0.89, 2.41, 6.03 x 10^° cm'^ respectively, all with very low velocity (-10 ±100 km s"^). 
From the emission line analysis they found a fiow velocity of —600 km s~^ and claimed that 
the X-ray absorption and emission originate in different regions. 

Ark 564 is also known to have a strong UV absorber, characterized by O vi , Si xil , 
Si IV and CiV absorption lines (j Crenshaw et al. 2002, Romano et al. 2002). 



3. Observations and Data Reduction 

Ark 564 was observed with the Chandra High Energy Transmission Grating Spectrom- 
eter (HETGS) on 2000 June for 50 ks and between 2008 August 26 to 2008 September 6 
for total duration of 250 ks. It was also observed with Low Energy Transmission Grat- 
ing Spectrometer (LETGS) on 2008 April 21 for 100 ks. Table 1 lists the observation log. 



Matsumoto et al.l analyzed and presented the results of the 50 ks Chandra HETG observa- 
tion done in 2000 June. Here we report on the Chandra archival 2008 observations of Ark 564 
made with both HETGS and LETGS. 

The HETGS consists of two grating assemblies, a high-energy grating (HEG) and a 
medium-energy grating (MEG). The HEG bandpass is 0.8 — 10 keV and the MEG bandpass 
is 0.5 — 10 keV but the effective area of both instruments falls off rapidly at either end of the 
bandpass. We performed the spectral analysis over 5 — 25 Arange. The LETG is combined 
with the Advanced CCD Imaging Spectrometer-Spectroscopic (ACIS-S) array or with the 
High Resolution Camera- Spectroscopic (HRC-S) array. The Ark 564 LETG observation was 
done with HRC-S array. The LETG/HRC-S has a band pass of 0.07-7.3 keV or 1.7- 170A, 
but due to the low S/N of data on either ends, we restricted our spectral fitting to 5 — 40 A. 

We reduced the data using the standard Chandra Interactive Analysis of Observations 
(CIAO) software (v4.3) and Chandra Calibration Database (CALDB, v4.4.2) and followed 
the standard Chandra data reduction thread^. For the Chandra ACIS/HETG observations. 



^http : //cxc ■ harvard . edu/ciao/threads/ index . html 
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we co-added the negative and positive first-order spectra and built the effective area files 
(ARFs) for each observation using the fullgarf CIAO script. Unlike ACTS, the HRC does not 
have the energy resolution to sort individual orders, and each spectrum contains contributions 
from all the diffraction orders. For the HRC/LETG observation we used the standard ARE 
files of orders 1 through 6 and convolved them with the relevant standard redistribution 
matrix files (RMF). 

For HETG observations, we generated the light curves in 2 ks bins for the energy band 
(0.3 — 10 keV), as shown in figure 1. The time-average count rate varies from 0.24 cts s~^ 
to 0.32 cts s~^ among 2008 observations. However, except for the initial 20 ks of obsID 
10575, the count rate of all observations are consistent with each other, with average value 
of 0.31 cts s-\ 

We analyzed the spectra using the CIAO fitting package Sherpa. As noticed above, the 
HETG observations do not show large variations, so to increase the signal to noise (S/N) 
of the spectrum we co-added the spectra obtained with each observation and averaged the 
associated ARFs using the ciao script add_grating_spectra. This gave a total net exposure 
time of 250 ks for the MEG and HEG. We fit the MEG and HEG data simultaneously, 
and discuss the LETG spectral analysis separately in section 5. Throughout the paper we 
applied the minimization technique in the spectral analysis and the reported errors are 
of Icr significance for one interesting parameter. 



4. HETG Spectral Analysis 

4.1. Continuum Modeling 

To model the intrinsic continuum of t he sou rce, we first fitted a simple absorbed (Galac- 
tic Nh = 6.4 X 10^° cm~^; iDickey et al.l (jl990l )) power law with varying photon index and 
amplitude. A single absorbed power law could not fit the data over the entire range. We 
found an excess of flux in the sp ectrum at energies below ~ 1 .5 keV. Ark 564 is known 



to have a strong soft continuum (ILeighly et al.lll999l . [Turner et al.lll999l ). so to fit this soft 
excess we added a black-body component to the above mentioned simple power-law. The fit 
improved significantly {x^/d.o.f. = 5882/3997, Ax^ = 515) and figure 2 shows the contin- 
uum model fit to the MEG spectrum. We also plot the data:fit residuals, which show strong 
absorption features consistent with the known WA of this source. 
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4.2. Local z = Absorption 



The spectra of Ark 564 show narrow abs orption hnes at ze ro redshift. The Ol and 
O II absorptions are attributed to the ISM ( ) Wilms et all boool ) and O vil , O vill and 
Fexvil absorptio n lines arise i n the circumgalactic medium (CGM) of our Galaxy or in 
the Local Group (iGupta et al.ll2012l . and references therein). We modeled all the statisti- 
cally significant local absorption features with Gaussian components of fixed width of 1 mA; 
the line response function of the grating is folded in the RMF file in each case. The model of 
continuum plus local absorption features ("Model A") improves the fit by Ax^ = 108 for 10 
fewer degrees of freedom. The measured equivalent width (EW) and statistical significance 
of each line are reported in Table 2 and labeled in figure 3. 



4.3. Intrinsic Absorption 

As reported above, the Ark 564 spectrum shows many strong intrinsic absorption fea- 
tures. We measured the position, EW and statistical significance of all intrinsic absorption 
lines by fitting negative Gaussians of fixed width of 1 mA(Fig. 3, Table 2). All the Ark 564 
absorption features are blueshifted with respect to the source, implying moderate outflow 
velocities of 82 — 239 km s~^. We find a few lines with no identification at 19.811 ± 0.006 A, 
19.850 ±0.005 A and 20.255 ±0.005 A in the observer frame. We marked these features with 
green underlines in figure 3. 



4.4. Photoionization model fitting: PHASE 



We used t he Photoionization mo del fitting code PHotoionized Absorption Spectral En- 
gine (PHASE; Krongold et al.l l2003l ) . to model the warm absorber features. The PHASE 
code self consistently models all the absorption features observed in the X-ray spectra of 
AGNs. At its simplest, an absorption-line spectrum can be fit with PHASE using only four 
input parameters: 1) the ionization parameter of the absorber U ; 2) the equivalent hydro- 
gen column density Nh; 3) the outflow velocity of the absorbing material Vout', and 4) the 
micro-t urbulent velocity Vty rh of the material. The abundances have been set at the Solar 
values ([Grevesse et al.lll993l ). We used the Ark 564 spectral energy distribution (SED) from 
Romano et al.l (120041 ) to calculate the ionization balance of the absorbing gas in PHASE. 
The SED constructed by Romano et al. is based on a quasi- simultaneous multiwavelength 
campaign, and thus is the most accurate overall SED for this source obtained to date. In 
the X-rays, however, we use our own fits, as this radiation is the one responsible for the ion- 
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ization of the charge states producing absorption in the X-ray band during our observations. 
We further connect the UV and the X-ray data with a simple power law (a straight line in 
log-log space connecting the last UV point and the first X-ray point). 

The most prominent absorption lines in the HETGS Ark 564 spectrum are those of 
Mgxi , Nex , Neix , Fexvil , Fexvill , Ovill , Ovil and Ovi at an outflow velocity 
of ~ —100 km s~^. We add a single ionized absorbing component to Model A described 
above (we call it Model B) to characterize this WA component. This absorber has best fit 
parameters of logU = 0.39±0.03, logNu = 20.94±0.02 cm~^, and an outflow velocity relative 
to systemic of —94 ±13 km s~^. The fit gives a significant improvement over model A (Fig. 4; 
X^/d.o.f. = 4722/3982, Ax^ = 1073). This absorber fits the high ionization lines so we will 
refer to this absorber as the "high- ionization phase (HIP)" component. The HIP component 
fits the absorption features produced by ion s such as Mg xi , Ne x , Ne ix , Fe xvill , Fe xvil , 



O VIII , and O Vll . iMatsumoto et al.l (120041 ) also report an absorber of similar characteristics: 
ionization parameter log f/ ~ 1 and absorption column Nh = 10^^ cm~^, derived using the 
column densities of O vill , Ne ix , Ne x and Mg xi . 

The single-absorber model does not fit the O vi absorption line and it also under-predicts 
the O VII absorption (Fig. 5, red curve). This suggests the presence of another absorber with 
lower ionization state. To fit the residual features we added another absorber to our previous 
model defining Model C. An absorber with ionization parameter logU = —0.99 ± 0.13, 
logNu = 20.11 ± 0.06 cm~^, and an outflow velocity of ~ — 137± 37 km s~^ successfully fits 
the low ionization lines, including O VI and O vil (Fig. 5 & Fig. 6). We call this component 
"lower-ionization phase (LIP)" absorber. This model gives a x^/d-O-f- = 4674/3979; Ax^ = 
48, significantly improving over the single-absorber model B. An F-test gives a higher than 
99.999% confidence for the presence of this absorber. The ionic column densities predicted 
by our best fit two-ionized absorber model (Model C) are listed in Table 4. As can be 
inferred from Table 4, except for O vil , the LIP and HIP components predict absorption 
from different ions but at similar velocities. This suggest that the Ark 564 absorbers may be 
present in the same outflowing multi-phase medium. We come back to this in section 6.3. 

To extrapolate the two component ionized absorber model of the MEG spectrum of 
Ark 564 to lower wavelengths, we simultaneously fit the HEG and MEG spectra (Table 3, 
Fig. 7). The HEG intrinsic absorption features are well fitted with HIP absorber, including 
Si XIII and Mgix lines which were not detected in the MEG spectrum due to low signal to 
noise. The best-fit model parameters of the MEG+HEG fit are consistent with the MEG- 
only fit. 



Though most of the Ark 564 intrinsic absorption features are well fitted with two warm 
absorbers, the model does not fit the unidentified absorption lines mentioned in section 4.3. 
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The identification, evaluation and interpretation of tliese features is be discussed in detail in 
a companion paper (Gupta et al. 2013b). 



5. LETGS Spectral Analysis 



Ark 564 was also observed with LETGS in April 2008. iRamirez fc Komossa I (120101 ) 
analyzed that observation and found a very weak absorption line of O vil Ka but a strong 
feature at 18.62 A (at the wavelength of Ovil K/3). For this reason they identified the 
absorption feature at 18.62 A as blueshifted Ovill Ka with velocity ~ 5500 km s~^. They 
further modeled the 17 — 25 A spectral region with two high ionization absorbing components 
with log{^) ~ 3, one at f ~ km s~^ and one at f ~ 5500 km s~^. Since we found no evidence 
of an outflow with velocity ~ 5500 km s~^ in the 2008 HETGS spectra, we reanalyzed the 
LETGS observation to check for the consistency with our model derived from the HETG 
data. 

To fit the Ark 564 LETGS spectrum continuum, we used the same model as for HETG 
data (a power law plus black-body). The best fit power law photon- index and temperature 
of black-body are reported in Table 3. We observed that the flux (2 — 10 keV) of the source 
varied from 2.48 ±0.11 x 10~^^ erg s~^ cm"^ to 2.79 ±0.15 x 10"^^ erg s~^ cm"^ between 
the HETG and LETG observations. To fit the intrinsic absorption features, we start with 
the two component ionized absorber model obtained for HETG spectra. This model fits the 
absorption features of Fexvill , Fexvil , Nex , Neix reasonably well, but overestimates 
the Ovil ka line strength {x^/do.f. = 1426/1190; Fig. 8 & Fig. 9). Allowing the PHASE 
parameters to vary freely, the fit improved considerably {x^/do.f. = 1370/1184; Fig. 10). 
The best fit parameters of LETGS WA model are reported in Table 3. This model fits 
the narrow absorption due to Ovil Ka, but leaves the residuals at 19.1 A, corresponding 
to O VII K(3 (Fig. 9) . We a lso tried to fit the LETG spectrum with models suggested by 



Ramirez &: Komossa I (|2010[ ). but the fit was not good {x /do.f. = 1936/1190). Though 



this model well fits the line at 19.1 A(observed frame), it also predicts other absorption 



lines which are inconsistent with the data. As noted by iRamirez fc Komossa I ( 120101 ). the 
absorption line at 19. lA is too strong to be by Ovil KP and could be, in part, due to a 
transient high velocity outfiow component (Gupta et al. 2013b, companion paper). 

Between the LETG and HETG warm absorber models, the HIP component parameters 
are consistent within errors. However, the LETG LIP has lower ionization parameter and 
higher column in comparison to HETG LIP. We also note that the column densities of highly 
ionized ions (O Vlll , Fexvil , Fexvill , Neix , Nex , Mgxi , Mgxil , and Sixil ) are higher 
while for those of less ionized ions (O vil and O VI ) are smaller in the HETG observation 
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than in LETG (Table 4). Both the instruments have good response in the spectral region 
where these lines are detected, so the observed differences cannot be due to instrumental 
artifacts; what is observed is the real variability in the WA properties between the two 
observations from 2008 April (LETG) to August /September (HETG). 



Discussion 



6.1. The Connection between UV and X-ray Absorbers 

Ark 564 was observed with HST (ST IS) and FUSE during May- July 2000 and June 2001 
respectively. I Crenshaw et al.l (|2002[ ) and iRomano et al.l ( 120021 ) detected intrinsic absorption 
lines blueshifted by ~ —190 krn . s~^ and ~ —120 km s~^ respectively, similar to that of 
the X-ray WAs. ICrenshaw et al.l modeled the UV data with a single absorber of ionization 
parameters logll of ~ 0.033 and column density logNn of ~ 21.2 cm~^. The UV absorption 
model also predicted the colum n densities of O Vll and Ovill of < 2.2 x 10^^ cm~^ and 



< 1.1 X 10^^ cm ^ respectively ( iRomano et al.ll2002l ). The Ovil column density measured 



from the X-ray WA models (0.3 — 1.3) x 10^'' cm~^ is consistent with UV upper limits. The 
total HIP+LIP O VIII column density of (2.0—2.9) x 10^'' cm~^ is an order of magnitude higher 
than the UV estimates. However, the LIP O Vlll column density = (3.0 — 4.6) x 10^^ cm~^ 
is in agreement with UV models. The consistency between the X-ray LIP absorber outflow 
velocity, hydrogen column density, and ionic column densities with the UV absorber model 
suggests that both are the same absorber. The HIP absorber, on the other hand, is different 
from the UV absorber, as expected. 



Romano et al.l from FUSE observations of Ark 564 measured the intrinsic O vi column 
density of (5.70 — 6.01) x 10^^ cm~^. The upper limit on Ovi column density measured 
from UV data is much smaller than our lower limit on Ovi column of 1.0 x 10^^ cm~^ in 
X-rays. This could be due to the variability of the WA between the two observations. We 
note, however, that X-ray and UV O v i column densities ha ve been found to be discrepant 
i n other AGN absorpt ion systems (e.g.. iKrongold et al.ll2003l ) and in redshift zero absorption 
Jwilhams et al.l 120061 ). 



6.2. Estimates on mass and energy outflows rates 

Using the values of warm-absorber parameters such as column density, ionization pa- 
rameter and outflowing velocities, we can give a rough estimate of the mass outflow rate 
[Mout] and the kinetic energy carried away by the warm absorbing winds (i.e. kinetic lumi- 
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nosity, Lk)- But before we can measure the mass and energy outflow rates, we must know 
tfie location of the absorber. However, in the equation for the photoionization parameter 
{U oc L/neR"^), the radius of the absorbing region (R) is degenerate with the density (ne). 
In principle we can put constr aints on the distance R by measuring the de nsity of the WA 



with variability analysis (e.g., iKrongold et al.l 120071 ). Using this technique, iKrongold et al. 



managed to determine the absorber density in NGC4051 and so the distance. For Ark 564, 
no such study is available in literature and in the present work the source does not show 
any significant variability either. Therefore, we will only set limits on the mass and energy 
outflow rates using the expression derived in IKrongold et al.l . Mout ~ i-'irrmpNHVout'r- 



The estimate of maximum distance from the central source can be derived assuming 
that the depth Ar of the absorber is much smaller than the radial distance of the absorber 
(Ar << r) and using the deflnition of ionization parameter {U = .'^j,^} J , i.e. r < Tr, 



2GMbh 

"ov.t 



a^unhc "'■^ several papers lower limit on the absorber distance was determined assuming that 
the observed outflow v elocity is large r than the escape velocity at r: i.e. r > 
However, as shown in iMathur et al.l (120091 ) WA outflow velocities are usually lower than 
the escape velocities, so cannot really be used to derive a lower limit on r. Using the best 
flt values of ionization parameter and column density, we estimated the upper limits on 
HIP and LIP absorber locations of thip < 40 pc and riip < 6 kpc respectively, which are 
not very interesting limits. Using the above equation and outflow velocities of 94 km s~^ 
and 137 km s~^, we obtain the mass outflow rates of Mout < 6.4 x 10^"^ g s~^ and Mout < 
2.2 X 10^^ g s~^ for HIP and LIP absorbers respectively. Similarly we obtained the constraints 
on kinetic luminosity of the outflows of Ek < 2.8 x 10^^ erg s~^ and Ek < 2.1 x 10'^° erg s~^ 



for the HIP and LIP absorbers respectively. 



^In comparison to the Ark 564 bolometric luminosity of 2.4 x 10^^ erg s ^ (IRomano et al. 

20021 ). the total kinetic luminosity of these outflows is Ex/Lhoi < 0.0001% for the HIP and 



< 0.006% for the LIP. Thus it is very unlikely that these outflows significantly affect the 
local environment of the host galaxy. The AGN feedback models typically required ~ 5% 
of the bolometric luminosity of an AGN to be conve rted into kinetic luminosity to have a 
signifi cant impact on the surrounding environment (jSilk &: ReesI Il998l IScannapieco &: Oh 



2004). 



6.3. Pressure Balance between LIP and HIP 



The presence of two different absorbing components with different temperatures but 
similar outflow velo city suggests that the absorber could arise from two phases of the 
same medium (e.g., lElvis et al.ll2000l IKrongold et al.l |2003[ ) . This is further supported by 



-li- 



the f act that multiple components of the ionized absorbe r are found in pressure balance 



Krongold et al.l l2003l: l2005l : 120071 ICardaci et al.ll2009l lAndrade- Velazquez et al.l 12010 



Zhang et al.ll2010[ ). T his result has proven valid against different methodologies and codes 
used in the analysis ( iKrongold et al.ll2013l ). To investigate whether or not the absorbing 
components of Ark 564 are in pressure balance , we generated the pressure-temperature equi- 
librium curve (also known as the "S-curve" ( iKrolik et al.lll98ll ). for the SED used in our 
analysis: Fig. 11). Interestingly, we find that the equilibrium curve of Ark 564 does not 
have the typical "S" shape where multiple phases can exist in pressure equilibrium, because 
there are no regions of instability (for all points in the plane the derivative of the curve is 
positive). 

Accepted at face value, this result implies that the absorber in Ark 564 is not in pressure 
balance and thus is not forming a multiphase medium, a result at odds with previous evidence 
on warm absorbers. While we cannot rule out this possibility, there are several arguments 
pointing towards a multiphase medium. We note that the two different absorbing components 
in Ark 564 have the same kinematics, which suggests that they are related. If they share 
the location (the most reasonable assumption), there must be a gradient of pressure between 
them, given that the LIP pressure is over 5 times larger than that of the HIP. Therefore, 
these two components should move on the S-curve to form a single component in a time 
comparable to the free expansion time, given by texp = Ai?/V^ (where AR is the thickness of 
the absorber and Vs the speed of sound in the medium). The flow time of the components (i.e. 
the time in which the components cross our line of sight to the source) is given to first order as 
tfiow = R/Vout (where R is the distance from the illuminated face of the absorbing material to 
the ionizing source and Vout its outflow velocity). For the warm absorber in Ark 564 Vs ~ Vout 
(specially for the HIP that is hotter). It follows that texp/ 1 flow = AR/R < 1. Then, the free 
expansion time is smaller than the flow time, and the two phase should dissolve into a single 
component before moving out from our line of sight, which is clearly not consistent with the 
data. Thus, if the two phases are not in pressure balance they should not be connected, and 
the similar kinematics in this, and in many other sources, would have to be considered a 
coincidence. 

Alternatively, the sources might be in pressure balance forming a multiphase medium, 
but there might be additional heating and/or cooling processes acting on the gas, changing 
the shape of the S-curve, but not the ionization b a lance . The most obvious parameter 
for this is the gas metallicity. iKomossa fc Mathun (120011 ) showed that the shape of the 
equilibrium curve not only depends upon the SED of the source, but also on the metallicity 
of the absorber, which affects the cooling of the gas. They further show that super-solar 
abundances restore the equilibrium zone in steep sp ectrum sourc e s and increase the pressure 
range where a multiphase equilibrium is possible. [Fields et al.l (|2007| ) showed that this is 
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indeed the case for the ionized absorber in Mkn 279. 



Since Ark 564 also has a steep spectrum and a monotonically rising "S" curve, we 
generated a new pressure-temperature curve with super-solar metallicity of 10 solar, shown 
as a dashed curve in figure 11 . This is a good as sumption as supersolar metallicity has 
been suggested for this source (IRomano et al.ll2004j ). Even with super-solar metallicity, the 
"S-curve" is very steep, without regions where multiple components can coexist in pressure 
equilibrium. This is due to fact that the continuum X-ray spectrum of Ark 564 is not only 
steep (r > 2.4), it also has an additional prominent soft exces s, similar to the behavior 
seen in sources with steep soft X-ray spectra (e.g., NGC 4051, iKomossa &: Mathurl 12001 



Mrk279, iFields et al.l 120071 ). In fact, additional modeling shows that the main reason for a 
steep "S-curve" is the extra heating produced by the soft excess (particularly in the LIP). If 
the soft excess continuum is not impinging directly on the absorbing gas, perhaps because 
it is the result of reflection toward our line of sight, then the two components would be 
in pressure balance. Other possibilities to produce a multi- valued S-curve, and LIP/HIP 
components in pressure balance include a weaker IR radiation field illuminating the material 
than the one observed (Krolik and Kriss 2001) or a dditional (more exo tic) sources of heating 
at high temperatures, such as those discussed in iKrolik et al.l ( 1l98ll ). Therefore, if warm 
absorbers are indeed a multiphase medium in pressure equilibrium, it is likely that the overall 
radiation field impinging on the gas is different than the one observed. This effect might be 
stronger in sources with steep soft X-ray spectra. We note, however, that this suggestion is 
speculative; we cannot prove it or rule it out. We also note that photoionization models of 
the Broad Line Region dera and that the ionizing continuum is different than the one observed 
(IBinette fc Krongoldll2008l . and references therein) . If warm absorbers are indeed in pressure 
balance, their S-curves can be used for a better understanding of the physical properties and 



the processes acting on the material ( iKomossa fc Mathurl l200ll . IChakravorty et al.l 12012 
Kroneold etliDboisi l. 



7. Summary 

Our best fit model of intrinsic absorption of NLSl galaxy Ark 564 requires a two- 
phase warm absorber with two different ionization states (HIP and LIP). Both the absorbers 
are outflowing at low velocities of order of ~ 100 km s~^. The HIP absorber reproduces 
most of the spectral features observed in the HETG spectra (Ovill , Neix , Nex , Mgxi , 
Fe XVII and Fe xvill ) except for a few at lower energies (O vil and O vi ), which are modeled 
by the LIP component. The pressure-temperature equilibrium curve for the Ark 564 warm 
absorber does not have the typical "S" shape, even if the metallicity is super-solar; as a 



- 13 - 



result the two WA phases do not appear to be in pressure balance. We speculate that the 
continuum incident on the absorbing gas is perhaps different from the observed continuum. 
We observe clear variability in the WA properties between the 2008 HETG observations 
and previous observations which could be in response to the change in continuum or the 
absorbing clouds passing our sight-line; the large time gap between observations does not 
allow us to distinguish between the two possibilities. 

We also estimate d the mass outflow r ate and associated kinetic energy assuming a 



biconical wind model ( iKrongold et al.l 120071 ) and flnd that it represents a tiny fraction of 
the bolometric luminosity of Ark 564. Thus it is highly unlikely that these outflows provide 
signiflcant feedback required by the galaxy formation models. 
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Table 1. Ark564 Chandra Observation Log. 



ID 


Start Time 


Exposure (sec) 


HETGS-ACIS-S 


9899 

9898 
10575 


2008-08-28 12:51:50 

2008-09-06 02:14:14 
2008-09-07 22:05:20 


84077 

99528 
62216 


LETGS-HRC-S 


9151 


2008-04-24X06:04:06 


99962 
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Table 2. Absorption lines observed in the Ark 564 Chandra HETG-MEG spectra. 





EW 




Ion Name ^ 


^rest 


Vout 


A 


mA 






A 


km s^^ 


9.393 ± 0.002 


7.1 ± 1.1 


44 


Mgxi 


9.17 


126 ± 64 


11.828 ±0.004 


5.5 ± 3.3 


25 


Ne IX 


11.54 


128 ± 101 


12.431 ± 0.001 


13.6 ± 1.4 


180 


Nex 


12.13 


82 ±24 


13.625 ±0.011 


5.0 ±1.8 


2 








13.775 ± 0.001 


15.2 ±0.8 


168 


Neix 


13.45 


107 ± 22 


14.165 ±0.005 


9.1 ±2.4 


22 


Fe xvii 


13.83 


126 ± 106 


14.554 ±0.004 


11.1 ± 1.6 


54 


Fe xviii 


14.21 


177 ± 82 


15.005 ± 0.005 


6.6 ±2.0 


10 


Fc XVII 


15.01 


local 


15.378 ± 0.003 


17.8 ± 1.0 


10 


Fe XVII 


15.01 


179 ± 58 


18.960 ± 0.004 


9.4 ±4.1 


8 


O VIII 


18.97 


local 


16.398 ± 0.002 


13.4 ±2.4 


70 


Oviii 


16.01 


91 ±37 


19.085 ± 0.005 


15.2 ± 2.6 


31 


Ovii /3 


18.63 


135 ± 79 


19.434 ± 0.004 


19.1 ± 3.2 


122 


Oviii a 


18.97 


111 ±62 


19.811 ± 0.006 


14.7 ±2.5 


27 








19.850 ± 0.005 


16.4 ±2.5 


34 








21.607 ± 0.007 


12.0 ±1.9 


27 


Ovii 


21.6 


local 


22.136 ± 0.008 


48.7 ± 7.0 


49 


Ovii 


21.6 


99 ± 108 


22.565 ± 0.005 


26.2 ±5.5 


24 


Ovi 


22.03 


239 ± 66 


20.255 ± 0.005 


14.8 ±4.0 


20 








23.288 ± 0.002 


19.1 ±4.7 


23 


On 


23.27 


local 


23.354 ±0.011 


25.4 ±5.8 


20 








23.506 ± 0.001 


22.5 ± 3.6 


40 


Oi 


23.45 


local 



^Change in after inclusion of Gaussian over continuum model 
''The " " marked the features with no identification 
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Table 3. Model parameters for the Ark 564 Chandra HETG and LETG spectra 



Units MEG MEG+HEG LETG 



Powerlaw 

Photon Index (r) 2.34 ± 0.04 2.38 ± 0.06 2.85 ± 0.05 

Normalization 10'^ ph keV'^ s'^ cm-"^ 13.4 ±0.3 13.7 ±0.5 20.8 ± 0.9 

Black Body 

kT keV 0.138 ±0.001 0.131 ±0.006 0.147 ±0.006 

Normalization^' 10"^ Lsg/Dfg 47.7 ± 1.6 49.3 ± 2.4 32.5 ± 2.3 

Warm Absorber: HIP 

Log U 0.39 ± 0.03 0.39 ± 0.04 0.48 ± 0.09 

Log N„ cm.-^ 20.94 ± 0.02 20.89 ± 0.03 20.78 ± 0.06 

Vout km s-i 94 ± 13 78 ± 15 83 ± 26 

Vturb km s-i 120 ± 150 160 ± 150 87 ± 150 

LogT^ K 5.20 ±0.01 5.20 ± 0.01 5.23 ± 0.07 

LogT/U{(x P f K 4.81 ± 0.04 4.81 ± 0.04 4.75 ±0.19 
Warm Absorber: LIP 

LogU -0.99 ±0.13 -1.04 ±0.22 -1.33 ±0.09 

Log Nh cm-^ 20.11 ±0.06 20.00 ±0.17 20.51 ±0.02 

Vout km s-i 137 ± 37 144 ± 64 189 ± 22 

Vturb km s-i 193 ± 150 186 ± 150 80 ± 150 

Log K 4.49 ± 0.02 4.48 ± 0.02 4.40 ± 0.09 

Log T/U {(X P)" 5.48 ±0.14 5.52 ±0.23 5.73 ±0.13 

4674 5131 1370 

Degrees of freedom 3979 6792 1184 

'^where L^^ is the source luminosity in units of 10^^ erg and Dio is the distance to the source in units 
of 10 kpc. 

''Derived from the column density and ionization parameter, assuming photoionization equilibrium. 

'^The pressure P oc rieT. Assuming that both phases lie at the same distance from the central source 
Tie oc 1/U, and P oc T/U. 
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Table 4. Ionic Column Densities Predicted By Models. 



Ion 




HETG 






LETG 




IoqNhip 


logNLip 


logNsuM 


logNnip 


IoqNlip 


logNsuM 


Caxiv 


14.86 




14.86 


14.67 




14.67 


Sxiv 


15.16 




15.16 


15.15 




15.15 


Caxv 


14.69 




14.69 


14.61 




14.61 


N VII 


16.17 


15.41 


16.24 


15.92 


15.36 


16.03 


Oiv 


10.34 


15.05 


15.05 




16.26 


16.26 


Ov 


12.46 


15.97 


15.97 


11.96 


16.80 


16.80 


Ovi 


14.09 


16.33 


16.33 


13.68 


16.88 


16.88 


Ovii 


16.44 


16.75 


16.92 


16.12 


16.88 


16.95 


O VIII 


17.37 


15.80 


17.38 


17.14 


15.53 


17.15 


Fe XVII 


16.06 




16.06 


15.90 




15.90 


Fe XVIII 


15.86 




15.86 


15.81 




15.81 


Ne IX 


16.61 


14.96 


16.62 


16.36 


14.37 


16.36 


Nex 


16.68 


13.18 


16.68 


16.54 


12.19 


16.54 


Mgx 


15.59 


13.63 


15.49 


15.21 


12.87 


15.21 


Mgxi 


16.37 


12.80 


16.37 


16.19 


11.77 


16.19 


Mgxii 


15.74 


10.34 


15.74 


15.67 




15.67 


Si XI 


15.75 


12.92 


15.75 


15.47 


10.80 


15.47 


Si XII 


15.79 


11.76 


15.79 


15.60 




15.60 


Si XIII 


16.14 




16.14 


16.05 




16.05 
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Fig. 1. — Ark 564 light curve of the HETG observations analyzed in this work, binned at 2 
ks resolution. Except for the initial 20 ks of obsID 10575, the count rates of all observations 
are consistent with each other, with average value of 0.31 cts s~^ (dash curve). 





5 10 15 20 25 




Fig. 2. — Top panel: The Ark 564 co-added Medium Energy Grating (MEG) spectrum in 
the observer frame. The red sohd hues show the best fit continuum model that consists of an 
absorbed power law and a black body component. Bottom panel: Plotted are the residuals 
showing strong WAs features. 
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Fig. 3. — Same as fig. 2; in addition to the continuum model, the absorption features are 
fitted with Gaussians. Note the numerous warm absorber features labeled in red, above the 
lines. The local (z=0) features are labeled in blue, below the lines. The unidentified features 
are marked with green underlines. The identification and interpretation of these lines are 
discussed in a companion paper. 
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Fig. 4. — Same as fig. 3, but the intrinsic absorption lines are modeled with PHASE. Only 
the high ionization phase (HIP) absorber is shown. 
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22 22.2 22.4 22.6 

Wavelength (X) 

Fig. 5. — An enlarged view of fig. 4 near tlie regions of O vil KP (top) and O vil Ka and 
Ovi (bottom). As can be observed, the Ovi line is not modeled by the HIP component 
(red curve) and it also underpredicts the O Vll absorption. The green curve shows the low 
ionization phase (LIP) of our model, reproducing the O VI and O Vll absorption. 



- 26 - 




Fig. 6. — Same as fig. 4, but showing both HIP (red) and LIP (green) components of Model 
C. The LIP only contributes at lower energies, particularly to O vil and O vi lines. 
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Wavelength (A) 



Fig. 7. — The High Energy Grating (HEG) spectrum of Ark 564 in the observer frame. The 
blue and red hues show the continuum and the WA model respectively. All the intrinsic WA 
features are well modeled with the HIP absorber. 
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Fig. 8. — The Low Energy Grating (LEG) spectrum of Ark 564 in the observer frame, fitted 
with Model-C of the HETG. Though most of the WA features are fitted well, but this model 
over estimates the Ovil Ka line at 22.13A. 
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Fig. 9. — An enlarged view of Fig. 8 near the region of O vil Ka (top) and O Vll KP (bot- 
tom). The red and magenta curve show the best fit HETG and LETG models respectively. 
The HETG model overestimates the O Vll ka (red curve, top panel), while the best fit LETG 
model underestimates O Vll k/3 (magenta curve, bottom panel) 
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Fig. 10.— 



Wavelength (A) 

Same as fig. 8, showing the LETG best fit two absorber model. 
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Fig. 11. — Pressure equilibrium curve (S-curve) for the Ark 564 SED used in the present 
analysis. The black curve is for the solar metallicity while the dashed blue curve is for super- 
solar (10 solar) metallicity. The points are for the LIP (lower) and HIP (upper) components. 



